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SUMMARY 


A method for assessing the wall interference in transonic two-dimensional 
wind tunnel tests has been developed and implemented in a computer program. 

The method involves three successive solutions of the transonic small distur- 
bance potential equation to define the wind tunnel flow, the perturbation 
attributable to the model, and the equivalent free air flow around the model. 
Required input includes pressure distributions on the model and along the top 
and bottom tunnel walls which are used as boundary conditions for the wind 
tunnel flow. The wall-induced perturbation field is determined as the differ- 
ence between the perturbation in the tunnel flow solution and the perturbation 
attributable to the model. The methodology used in the program is described 
and detailed descriptions of the computer program input and output are pre- 
sented. Input and output for a sample case are given in an appendix. 


INTRODUCTION 

A method has been developed for assessing the wall-induced velocity pertur- 
bation existing during transonic airfoil tests in two-dimensional wind tunnels. 
The method can be considered as a transonic nonlinear counterpart of the 
classical wall interference theory. The basic principles are outlined in refer- 
ence 1 and a qualitative description of the method is given in reference 2 along 
with a discussion and illustration of its capabilities. The present implementa- 
tion of the method includes some capabilities which were added subsequent to the 
date of reference 2. 

The present paper is intended to serve as a user's guide to the computer 
program TWINTAN* in which the method has been implemented. This paper presents 
the methodology used in the program and a detailed description of the input 
required and the output provided by the program in its present form. Input and 
output for a sample case are given in an appendix. 

The computer program described herein is in an elementary form suitable 
for evaluation and exploratory use. Current experience with the program leads 
to the belief that the Mach number correction is generally reliable and that 
the nonuniformities in the wall-induced perturbation field and their effect on 
the airfoil pressure distribution are adequately revealed. The accuracy of the 
angle of attack correction, however, depends directly on the accuracy with which 
a flow direction reference for the wind tunnel flow computation can be estab- 
lished. It is anticipated that improvements in both accuracy and convenience 
could result from mutual optimization of this program with the instrumentation 
and data system of the particular test facility with which it is to be used. 

^Program TWINTAN is distributed through COSMIC, Computer Software Manage- 
ment and Information Center, 112 Barrow Hall, University of Georgia, Athens, 

GA 30602. 
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SYMBOLS 


a coefficient of (py in boundary condition expressions 

A,B ,C,D coefficients in tridiagonal form of governing equation 

B 1 , D* coefficients in line relaxation algorithm 

c airfoil chord 

c section lift coefficient 

1 

f transformation derivative in longitudinal direction 

g transformation derivative in transverse direction 

j longitudinal grid index, also longitudinal computational coordinate 

k transverse grid index, also transverse computational coordinate 

M Mach number 

n stretching parameter in coordinate transformation 

P,Q quantities in the finite difference formulation of the governing 

equation 

q dynamic pressure 

s distance along airfoil surface 

U longitudinal velocity of uniform flow 

V resultant velocity 

x longitudinal physical coordinate 

y transverse physical coordinate 

z general physical coordinate 

a angle of attack 

r circulation 

Y ratio of specific heats 

A prefix denoting an increment 

6 scaling parameter in coordinate transformation 
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£ general computational coordinate 

A coefficient of longitudinal term in governing equation 

y switching parameter in finite difference formulation of the governing 

equation 

$ total velocity potential 

<j) dimensionless perturbation velocity potential 

Subscripts: 

x,y denote partial differentiation with respect to x or y 

j,k denote value at j-th or k-th grid line 

R wind tunnel reference conditions 

00 corrected far field conditions 


METHODOLOGY 

The procedure implemented in program TWINTAN involves three successive 
finite-difference solutions to the two-dimensional transonic small disturbance 
potential equation. The first solution is a calculation of the flow around the 
test airfoil in the wind tunnel. Pressure distributions measured during the 
test to be assessed, both on the airfoil and at the upper and lower walls, are 
imposed as boundary conditions to enhance the fidelity of reproduction of the 
actual tunnel flow. The assumptions and inaccuracies usually encountered in 
modeling the ventilated walls and the remaining tunnel geometry are thereby 
avoided. 

The second solution is a calculation of the flow perturbation attributable 
to the model. Free air outer boundary conditions are used and the model bound- 
ary conditions are imposed in the form of singularity distributions extracted 
from the tunnel flow solution. The wall-induced perturbation is then determined 
as the difference between the perturbations calculated in the two flow solutions. 
In order to keep the nonlinear effects comparable, the streamwise component of 
the wall-induced perturbation at a user specified match point is evaluated 
iteratively during the second solution and is transferred as a correction to the 
far field velocity. Upon convergence of the second solution, the wall-induced 
corrections to Mach number and free-stream static and dynamic pressures are 
thereby fixed. 

The third solution is a calculation of the free-air flow at the corrected 
free-stream Mach number around the equivalent airfoil shape determined in the 
tunnel flow solution. The test lift coefficient, corrected by the dynamic 
pressure ratio, is specified and the angle of attack adjusted to satisfy the 
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Kutta condition. The results include a specific value for the angle-of-attack 
correction and a pressure distribution which could be considered fully corrected 
for wall interference if it does not depart excessively from the original air- 
foil pressure distribution reduced according to the corrected free-stream con- 
ditions. Some features of the method are described in detail in the following 
sections to aid the user in adapting the program to a specific facility as well 
as in making refinements. 


Governing Equation 

The dependent variable in each of the flow solutions of the present program 
is the dimensionless perturbation potential (J) which is related to a total 
dimensional velocity potential $ by 

* = (1) 

where U R is the velocity corresponding to the tunnel reference Mach number. 

For the tunnel flow solution is set equal to U R but during the model 

perturbation solution is updated by adding the wall induced blockage 

velocity. The updated value is used for the equivalent free air solution. 

Thus, by equation (1) <j> is scaled to the same units in all solutions while 

representing the perturbation from uniform flows of different velocities. 

The perturbation potential is governed by an extended form of the transonic 
small disturbance equation expressed as 



(2) 


where 



(Y + 1) 


2 H* + ^r_ 

' 3x ^ 2U ro 3 x) 


The higher order term in A is included to improve the accuracy of velocity 
equivalence between the perturbation and the far field and thereby suppress 
spurious errors which were otherwise observed when U R ^U oo differed signifi- 
cantly from unity. 


The conservative form finite difference approximation to equation (2) used 
in the program is described by the following expressions: 
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In defining the transformation derivatives f and g, the quantities j and k 
each serve the dual role of computational variable and grid index because the 
transformation is scaled so that the mesh size is one unit of the computational 
variable. The first derivatives of $ appearing in the preceding expressions 
are all formed by single mesh differencing; for example: 


(1^,** * 9k+ ‘ I (' , ’j' k+1 " * j ' k ) 

Note that at shock points, ]ij = 0 and Pj-1 = 1 so the longitudinal term of 
the governing equation contains the sum of central difference and upwind differ- 
ence contributions. If the nonconservative form is selected (ICNSV = 0) the 
upwind difference contribution is omitted at shock points. 

In the following discussion, points in the lower half plane will be identi- 
fied by a negative sign on the k index. Accordingly, g_ k “ -g^- 
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The governing equation is arranged in the tridiagonal form 


\ ‘f’jjk-l + B k $j # k + c k ^k+l “ D k 


( 3 ) 


along j = constant lines and is solved by the method of successive line over- 
relaxation subject to boundary conditions imposed in Dirichlet form at the 
outer boundary and in various forms at the airfoil. The development of these 
boundary conditions for the three flow computations is described in subsequent 
sections. The implicit algorithm used to solve equation (3) along each line is 
based on the recursive relationship 


p j,k 


= D' 


+ B' 


k $j,k-l 


(4) 


where the coefficients B' and D' are defined by 


_ D k ~ C k D 'k+1 
° k " B k + c k B' k+1 

r| = “ A k 

B k " B k + c k B 'k+1 


and are calculated recursively starting with the outer boundary condition and 
proceeding inward to k = ±2. 

Because of the symmetric grid indexing used, the k = 1 grid line is 
treated numerically as a boundary along its entire length. At points ahead of 
the airfoil and in the wake, the boundary condition is derived from the 
governing equation. In the wake, the potential jump relation 


♦j.-l = *3,1 - r < 5 > 

is used along with equation (4) applied at k = ±2 to write equation (3) at 
k = 1 and obtain the solution. 

D 1 - C l [ p, 2 * °’-2 + F ( 1 - B '- 2 ?] 

* 3,1 B 1 + c l ( B ' 2 + B ' — 2 ) 


Equations (5) and (4) are then used to complete the line solution. The same 
procedure, but with T = 0, is used ahead of the airfoil. 
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Boundary Conditions for Tunnel Flow Computation 

All input pressure coefficients are interpolated onto the appropriate 
longitudinal mesh center points, converted to the form of velocity magnitude 
squared using the exact pressure coefficient definition, and stored in this form. 
The boundary conditions for the tunnel flow computation are derived primarily 
from these data. 

On the longitudinal outer boundaries at k = ±KW, values of <J> X are 
found from 



(7) 


and integrated to determine the boundary values of <f). Values of 4> y are 
assumed initially to be zero and then are updated iteratively. The 0 X gradi- 
ents approaching the ends of the upper and lower boundaries are used with the 
governing equation (2) to find <f>yy at each boundary corner. At the two 
upstream comers, <j> y is specified by input quantities. The variation of <J> 
on the upstream boundary is described as a fourth degree polynomial in y 
whose coefficients are determined from the values of <f)y and noted above 

and the assumption that cf) = 0 at the lower upstream corner. On the downstream 
boundary, <J> is defined as a cubic in y with the term T/2 sgny added to 
provide the needed potential jump across the wake. The four coefficients are 
determined to satisfy the known values of (j> and <j)yy at the comers. 

The airfoil boundary conditions are imposed in Dirichlet form on a slit 
at y = 0 rather than on the airfoil surface. It is not appropriate to 
require that the airfoil surface velocities be reproduced on the slit. Instead, 
we require that the longitudinal distribution of velocity potential on the 
airfoil surface be reproduced on the slit. Thus, 


/ 


ds 

V -r— dx = 
dx 


Iff + ♦*) 


dx 


( 8 ) 


where V is the airfoil surface velocity determined from the pressure coef- 
ficients and ds is the element of surface arc length. The present program 
calculates (f) x by equating the integrands of equation (8) using the approxi- 
mation 



in which <t>y is calculated at the slit and is updated iteratively. This 
procedure may be recognized as an adaptation of the Riegels rule. 

The approximations inherent in the above procedure are unacceptable at a 
blunt leading edge, therefore equation (8) is used to determine <|) x for only 
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those mesh intervals downstream of the leading edge. In the mesh interval 
bracketing the leading edge, upper and lower surface values of <j> x are 
determined such that two broad constraints are satisfied. The first requires 
that the total circulation determined by integrating A<J> X over the entire 
airfoil chord must correspond to the prescribed lift coefficient. The circu- 
lation error is corrected directly by adjusting the &(p x in the leading edge 
interval. The second constraint is the airfoil thickness closure condition 
which requires that the net source strength determined by integrating A<py 
over the entire airfoil chord must correspond to the prescribed drag coefficient. 
This constraint is implemented by allowing the error in net source strength to 
control the average value of <p x in the leading edge interval. The direct 
effect of negative <f> x at the leading edge is to depress the value of <J> at 
the first grid point behind the leading edge, thereby creating source strength 
at this point. As a secondary effect, however, added source strength is 
created over the entire chord with a distribution such that the (j) x boundary 
conditions are still satisfied in each mesh interval. The existence of such 
a source distribution which has no effect on <J> X anywhere on the chord line 
causes the thickness design problem to be ill posed without a specific con- 
straint on thickness closure. The role of the thickness closure constraint 
in the airfoil thickness design problem is analogous to the role of the Kutta 
condition in the airfoil lift analysis problem. 


Boundary Conditions for Model Perturbation Computation 

The purpose of the model perturbation computation is to define that part 
of the tunnel flow perturbation which is attributable directly to the airfoil. 
Accordingly, the airfoil boundary condition in this computation should assure 
that the vorticity and source distributions on the slit representing the air- 
foil duplicate those existing in the tunnel flow solution. Thus, the quan- 
tities 


A <j> . = <f> . 1 

J Drl 


- * 1.-1 


(9a) 


AcJ> . = <p 
yj 


- 

yj ' 1 yj /- 1 


(9b) 


as extracted from the tunnel flow solution are to be reproduced in the free 
air computation. 

The values of < py used in (9b) are those at grid points on the slit. 
They may be expressed in terms of the <p values at surrounding grid points 
by using a 4>yy at the slit expressed by the governing equation (2) . The 
results can be written in the following form: 


B 1 + C *),2 - D 1 ' a W (10a) 

B -1 *j,-l + C *y,-2 * D -1 + a (10b > 
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From equations (10) , (9) and (4) we find 

D + D - CU>' + D' ,) + (b + C B' 0 ^A(f) . -aA<j) . 

* , -1 =1 512 -2 L Lzl II (11, 

B i + B -i + c ( b ' 2 + b '- 2 ) 

which is used along with equation (9a) to impose the airfoil boundary conditions 
in the model perturbation computation. 

The values of the perturbation potential imposed on the outer boundary 
of this computation are calculated from a far field analytic expression which 
includes the two-dimensional lift and thickness integral terms given in 
reference 3 and an additional term representing the far field influence of the 
net source strength at the airfoil. 

During this computation, the far field Mach number is updated by trans- 
ferring the wall-induced longitudinal velocity perturbation evaluated at the 
specified match point into the far field velocity using the following procedure: 


U 

00 

— = 1 + 


u. 


R 



(12) 


where T is the value of 4> x at the match point in the converged tunnel 
flow solution and (4 >x)m that i* n current iteration step of the model 
perturbation solution. The updated Mach number is then found from 



The linearized compressibility factor used in 
conditions is updated along with M . 


(13) 


calculating the outer boundary 


Boundary Conditions for Equivalent Free Air Computation 

The equivalent free air computation is to define the free air flow at the 
corrected far field Mach number around the equivalent airfoil shape described 
by the tunnel flow solution at an angle of attack corresponding to the lift 
measured in the experiment. The lift coefficient based on the free stream 
conditions at M is 

oo 
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Ci = Cl — 

l °° lR ^co 


where is based on tunnel reference conditions and 


%> 


q 


R 




The far field boundary conditions and the potential jump in the wake are 
determined using 


U q 

CO 

r = r — — 
00 R 


where the velocity ratio accounts for the fact that T (like <f>) is expressed 
in units of U R c. 

The airfoil boundary condition is a Neumann condition using values of (J)y 
extracted directly from the tunnel flow solution by use of equation (10) . The 
boundary condition is then formed by combining equation (10) and (4) and 
allowing for an angle of attack correction A°c. 

D ±1 - C D '±2 ? a (Vi,±l " *“) 

+1 = 

“ B ±1 + C B ’±2 

where the upper and lower signs refer to the airfoil upper and lower surfaces, 
respectively. 


DESCRIPTION OF INPUT 

Program TWINTAN is written in FORTRAN 4 and has been implemented on the 
CDC CYBER 175 computer with a central memory requirement of approximately 
60,000 (octal) of 60 bit words. The program input is provided in two namelist 
blocks. The first, labeled NPUT, contains 12 parameters to define the compu- 
tational grid, 6 parameters for problem constraint arid 8 parameters for 
numerical process control. The second block, RDPX, provides arrays of wall 
and model pressure coefficients and their longitudinal locations. All input 
data are nondimensional with the airfoil chord assumed to be the unit of 
length. 
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Computational Grid Inputs 


The perturbation potential is stored by the program in an array addressed 
as PHI ( J , K, L) . The index L takes on values of 1 or 2 for the upper or lower 
half plane respectively. J and K are the longitudinal and lateral indices 
respectively in each half plane. The physical configuration of the upper half 
grid is sketched in figure 1. The K = 1 grid line is the axis of grid 
symmetry between the two half planes and represents the tunnel center line, 
y = 0. The lower half grid is a mirror image of the upper. The longitudinal 
grid spacing is uniform from Jl to J4. The longitudinal spacing upstream from 
J1 and downstream from J4 and the lateral spacing of the entire grid are 
stretched by transformations of the form. 



where z represents a physical coordinate and £ is the computational 
coordinate in the corresponding direction. In the present program, the trans- 
formation is truncated to prevent z reaching infinity. Equation (14) is used 
to relate the coordinate x or y directly to the grid index J or K by 
assuming a computational mesh interval A£ of unity. The stretching parameter 
n then represents the number of mesh intervals which would cover the range 
0 <_ z <_ °°, and the scaling parameter 6 is the limit as z approaches zero 
of the physical mesh interval Az. The specific relations for the upstream 
and downstream, and lateral transformations are formed by substituting the 
parameters tabulated below into equation (14) . 


Index Range 

1 1 j i J1 

J4 £ J < J5 

1 £ K £ KB 

K 

Jl - J 

J - J4 

K - 1 

z 

x(Jl) - X 

x - x ( J4) 

y 

n 

M 

M 

N 

5 

DXMN 

DXMN 

DYMN 


The program sets the uniform mesh interval Ax in the range from Jl to J4 
equal to DXMN, locates J3 at the airfoil trailing edge (x = 1.), and determines 
J2 such that the airfoil leading edge (x = 0.) lies in the mesh space between 
J2 and J2 + 1. The boundaries of the uniform longitudinal mesh region are 

x(Jl) = 1. - DXMN * (J3 - Jl) 

x ( J4) = 1. + DXMN * ( J4 - J3) 
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The foregoing information is sufficient to determine the physical coordi- 
nates x,y of each grid intersection J,K in terms of the following program 
input parameters. The longitudinal grid used in the free-air computation is 
defined by the integers Jl, J3, J4, J5 and M and the real constant DXMN 
which must satisfy the following constraints. 

J5 < 100 

Jl < M 

J5 - J4 < M 

{ J3 - Jl) * DXMN >_ 1. 

J4 > J3 

The grid lines serving as the upstream and downstream boundaries of the tunnel 
flow computation are denoted by the integers JU and JD respectively which 
must satisfy: 

JU £ 1 

JD <_ J5 

In addition, all wall boundary points, located at mesh centers from x(JU + 1/2) 
to x(JD - 1/2) must lie within the longitudinal range of the input wall 
pressures. 

The lateral grid is defined by the integers KB, N and KW and the real 
constant DYMN which must satisfy: 

KB £ 50 

KB £ N 

KW £ KB 

KW identifies the lateral boundary of the tunnel flow computation and must be 
selected along with N and DYMN such that ±y(KW) is the lateral distance from 
the tunnel center line to the upper or lower line of input wall pressures. 

The following points should be considered also in selecting grid input 
parameters. The outer boundaries of the free-air computation, x(l), x(J5) 
and y{KB) should lie several chord lengths from the model for best accuracy of 
the free-air boundary conditions. If, however, Jl and (J5-J4) closely approach M, 
and KB closely approaches N, the solution convergence rate is adversely 
affected. 
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Problem Constraint Inputs 


Five real and one integer input parameters provide constraints needed in 
the tunnel flow computation. AMREF is the reference Mach number corresponding 
to the tunnel reference conditions used in reducing the input data to coefficient 
form. The lift coefficient CL defines the total airfoil circulation. The 
accuracy of the calculated wall-induced upwash depends on the accuracy of the 
CL input. CD is the airfoil drag coefficient which provdes the airfoil 
thickness closure constraint. The integer JP defines the location 
(x = x(JP - 1/2) , y = 0) at which the local Mach number in the model perturba- 
tion computation is matched to that in the tunnel flow. For low Mach number 
cases, the results are essentially independent of the choice of JP. For 
supercritical flows, the match point should be located just upstream of the 
sonic point in the shock, thereby matching shock strength in the two flows. 
Numerical divergence might be encountered if the match point is just downstream 
of a strong shock. 

The real parameters SLA and SLB are the onset flow streamline slopes at the 
upper and lower corners respectively of the upstream boundary of the tunnel 
flow computation x(JU). The present formulation presumes that this upstream 
boundary is located in the solid nozzle region of the tunnel where the flow 
direction is constrained by the known wall slopes SLA and SLB. If this 
condition cannot be met or if other conditions, such as poor resolution of 
the wall pressures, prevent accurate modeling of this upstream region, SLA and 
SLB may be input as zero and the tunnel flow direction monitored by comparing 
the orientation of the airfoil shape calculated in the tunnel flow solution 
with the known incidence of the test airfoil. It has been observed that for 
supercritical cases, the calculated airfoil orientation is affected by the 
choice between conservative and nonconservative differencing at shock points. 
Because it is not yet clear which choice leads to the most accurate angle of 
attack correction, the choice is made optional as described in the following 
section. 


Computational Control Inputs 

The parameters discussed in this section affect the differencing form, 
convergence rate, stabilization, termination and amount of diagnostic output 
for the iterative process in each of the three solutions. Default values of 
each of these parameters are established in the program but may be overridden 
simply by including new values in the NPDT namelist. 

The integer parameter ICNSV, if left as its default value of 1 provides 
for fully conservative differencing of the governing equation. Nonconservative 
differencing at shock points may be selected by setting ICNSV = 0. 

The real parameter RFM affects only the model perturbation solution and 
is the underrelaxation factor used in transferring the Mach number correction 
evaluated at the match point to the far field. The default value of 0.2 is 
satisfactory for most cases but a smaller value might be needed to avoid 
oscillatory instability if the match point is located in or just upstream of a 
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very strong shock. The real parameter RFA, used only in the equivalent free 
air solution, similarly controls the angle of attack adjustment in response to 
a mismatch between the circulation integral around the airfoil and the 
circulation value prescribed for the far field and wake. 

The five remaining parameters are dimensioned quantities denoting values 
stored in three-element arrays. The array subscripts from 1 to 3 correspond to 
the three successive flow solutions. The real array ORF contains the over- 
relaxation factor used at subsonic points in the successive line overrelaxation 
procedure. The real array DXT contains a damping factor to retard updating of 
longitudinal velocity perturbations. The real array EPS contains the convergence 
criterion in terms of the absolute value of the perturbation potential change 
in one iteration at the grid location where this change is greatest. Iteration 
is termined either when this criterion is met or when the number of iteration 
steps exceeds the value obtained from the integer array IMX. information 
describing the convergence history is available for output at each step of 
iteration. The values stored in the integer array ITRC determine how much of 
this information is actually printed. An ITRC value of 1 causes printing at 
each iteration step; a value of 2 causes printing every tenth step and at 
termination; and the dafault value of zero causes only the final values to be 
printed so that it may be determined whether the iteration converged or ran 
the maximum number of steps. If either of the first two solutions is terminated 
by reaching the maximum number of iterations and the final value of the con- 
vergence parameter is greater than 10 times the criterion given in EPS, then 
all further processing of that case is aborted. 


Pressure Distribution Inputs 

In the tunnel flow computation, the outer boundary conditions are estab- 
lished from the input pressure coefficient distributions along or near the 
upper and lower tunnel walls. Similarly, the inner boundary conditions are 
established from input pressure coefficient distributions on the airfoil upper 
and lower surface. These data are input through parameters in the namelist 
block RDPX. The wall pressure data are described in a two-element integer 
array and two real two-dimensional arrays. The integers NW(1) and NW(2) are 
the numbers of pressure input locations along the upper and lower walls 

respectively- The first real array XW inputs the x coordinates of the wall 

pressure locations in order from upstream to downstream. The first index 
numbers these locations from 1 to the appropriate NW element and the second 
index is 1 or 2 for the upper or lower wall respectively. The second real 
array CPW inputs the pressure coefficients with the same order and indexing 
as their corresponding locations. The airfoil upper and lower surface pressure 
coefficients and their chordwise locations are input in a fully analogous 
manner using the integers NM(1) , and NM(2) and the arrays XM and CPM. All 
four real arrays are dimensioned to accommodate values of NW or NM up to 30. 

All input pressure distributions are interpolated onto boundary points 
located at the centers of the longitudinal mesh intervals. These boundary 

points range from J2 + 3/2 to J3 - 1/2 on the airfoil (K = 1) and from 

jU + 1/2 to JD - 1/2 on the outer boundary (K = KW) . The interpolation routine 
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fits a spline to each distribution with a zero second derivative condition at 
the spline ends. Because the routine as presently implemented does not provide 
for extrapolation, each input pressure distribution must cover the full range 
of boundary points. Manual extrapolation of the experimental pressures at 
the airfoil trailing edge might be necessary in some cases to provide the 
required range of data. Near the leading edge, on the other hand, the number 
of experimental measurements might far exceed the requirements of this program 
and, therefore, one might omit some of the experimental data. Care should be 
taken not to introduce spurious oscillations in the spline fit between data 
points. One might even input a false pressure coefficient at the leading edge 
to help the spline provide a smooth fairing in the range of interest. 


DESCRIPTION OF OUTPUT 

Output is generated during each major phase of computation so that the 
user may scrutinize the processes leading to the final wall-induced perturbation 
results. The process output will be described in four sequential phases, input 
processing, tunnel flow solution, model perturbation solution, and equivalent 
free air solution. All output velocities are expressed in units of the 
reference stream velocity corresponding to the input tunnel Mach number AMREF. 
All output velocity potentials are expressed in units of the product of this 
reference velocity and the airfoil chord. 


Input Processing 

The contents of the namelist block NPUT are listed first. The longitudinal 
grid structure is then described by listing the x location and the derivative 
dJ/dx. These quantities evaluated at the grid lines J are identified as 
X ( J) and Fl ( J) . The corresponding quantities evaluated at J + 1/2 are given 
by XMID(J) and F2 (j) . The lateral grid structure is described by the grid line 
location Y (K) , and the derivative dK/dy evaluated at the grid lines G1 (K) and 
at the K + 1/2 points G2 (K) . The user should check to see that Y (KW) is the 
lateral location of the input wall pressures. The results of interpolating 
the airfoil pressure coefficients onto the mesh midpoints XMID(J) are then 
listed in separate groups for the upper and lower surfaces. This listing can 
be used in conjunction with the input data to judge the quality of the spline 
fit. 


Tunnel Flow Solution 

The convergence history of the tunnel flow solution is shown by listing 
the iteration number I, the largest change in the perturbation potential from 
the previous iteration DPHIMX / and the grid intersection JMX, KMX where this 
change is located. Negative values of KMX indicate points in the lower half 
plane. The flag ISUP is 1 if this grid point is a supersonic point and 0 if it 
is subsonic. The number of supersonic points NS and the current value of air- 
foil trailing edge thickness TTE are also listed. Both the inner and the 
outer boundary conditions depend on the solution and therefore are updated 
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iteratively. To avoid using poorly developed data for the update, it is 
performed only when the convergence parameter DPHIMX is less than some threshold 
value. The threshold levels specified for the inner and outer boundary condi- 
tion updates are 10 ”^ and EPS (1) *10. respectively. The threshold levels are 
separated so that numerical problems associated with each update process might 
be identified more readily. Because the initial approximation to the inner 
boundary conditions is poor, it is normal for the convergence history to show 
that DPHIMX bounces up from the 10“3 level several times before sinking through 
it. This listing will end either with DPHIMX £ EPS(l) or with I - IMX(l) 
depending on whether or not convergence was achieved. 

The data in the following listings are identified by a heading as 
pertinent to the tunnel flow. The critical pressure coefficient is given as 
CPSTAR. Data evaluated at the airfoil surface are listed next. Two lines are 
printed for each grid index J. Data on the first line pertain to the uppper 
surface or both surfaces and those on the second line to the lower surface 
or the increment between surfaces. The first four columns give data evaluated 
at the center of the mesh interval upstream from J. The calculated shape of 
the airfoil surface streamlines is given by the coordinates X(J-.5) and YMOD. 

The pressure coefficient CP (REF) can be compared with the previous listing of 
interpolated pressure coefficients to see how well the inner boundary condi- 
tions were satisfied. The longitudinal perturbation velocity component U is 
obtained directly from the potential gradient along the K = 1 grid line. The 
data in the last four columns are evaluated at the grid points located at X(J) . 

V is the lateral component of perturbation velocity and DELTA PHI and DELTA V 
are the increments between the upper and lower surface values of the perturba- 
tion potential and the lateral velocity respectively. These increments are 
stored for use as the model boundary conditions in the model perturbation 
computation . 

The perturbation velocity components along the outer boundaries of the 
tunnel flow computational region are given next. The next listing pertains 
to the upstream and downstream boundaries. The different coordinate locations 
given for the two components reflect the fact that U is evaluated over longi- 
tudinal mesh intervals and V over lateral mesh intervals. In the case of the 
upper and lower boundaries listed next, the V component is extrapolated to grid 
points on the K = KW boundary. Thus, the listed values of U and V are given at 
the indicated longitudinal location along the upper and lower boundaries at 
K = KW. 

The final listing from the tunnel flow solution gives the perturbation 
potential at grid points in the immediate vicinity of the airfoil. The grid 
points are those on the K = 1, 2, and 3 grid lines in the upper and lower half 
planes for J from Jl to J4. 


Model Perturbation Solution 

Data output listed for the model perturbation solution is structured 
similarly to that described for the tunnel flow solution. This discussion, 
therefore, will highlight the differences between the two. 
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In the model perturbation solution, the far field Mach number is updated 
after each iteration cycle for which the convergence parameter DPHIMX is 
less than EPS (2) *10. As was previously indicated, the numerical stability 
of this update can be affected by the input parameters JP and RFM. The update 
stability characteristics can be judged from the convergence history by 
comparing the sign and magnitude of the UFAC changes immediately following 
two successive penetrations of the EPS (2) *10. threshold. 

After convergence, the final far field Mach number is given as AMINF, 
the corresponding velocity as UINF/UREF and the dynamic pressure ratio as QFAC. 
This far field condition represents the corrected reference condition for the 
tunnel data. Accordingly, the model surface pressures listed as CP (INF) as 
well as the critical pressure CPSTAR are in the form of pressure coefficients 
based on the corrected reference. The change from CP (REF) in the tunnel flow 
to CP (INF) at the match point JP - 1/2 is due simply to this reference change 
but, at other locations, it is also influenced by nonuniformity of the wall- 
induced blockage velocity. Similarly, the airfoil shape given by YMOD, 
when compared with that in the tunnel flow, might show not only a different 
angle of attack due to wall-induced downwash, but also differences in mean 
camber shape due to nonuniformity of wall-induced downwash. It should be noted 
that the perturbation velocities and the perturbation potential in the free-air 
solution output are perturbations from a uniform flow at the far field velocity 
UINF. The velocity unit, however, is still the tunnel reference velocity UREF. 

The wall-induced velocities are calculated as the difference between the 
perturbation velocities calculated at the same locations in the two flow 
solutions. This is performed in the present program only along the tunnel 
center line (K = 1) from JU to JD. The lateral velocities in the free air 
solution VF and due to the walls VW are listed at the grid point locations XV. 
The longitudinal velocity perturbations in the free air solution UF and due 
to the walls UW are given at the mesh center location XU. Over the model 
chord, values corresponding to both upper and lower surfaces are given. 

Although the wall-induced velocities should be identical on the two surfaces, 
some differences can appear in VW which reflect the accuracy with which the 
velocity jump boundary conditions are satisfied. 


Equivalent Free Air Solution 

The convergence history listing for the equivalent free air solution 
includes the current value of DA, the angle of attack correction in radians. 
This solution should converge in a relatively small number of iterations. If 
convergence difficulties are encountered, the iteration process can generally 
be stabilized by some combination of reduced RFA, reduced ORF(3), and increased 
DXT ( 3 ) . 

After convergence, the final value of the angle of attack correction is 
listed and results at the airfoil surface are given in the same form as those 
from the previous solutions. The airfoil shape described by YMOD should be the 
airfoil shape from the tunnel flow solution rotated through the angle of attack 
correction. The numerical procedure used, however, applies the angle of attack 
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correction as a surface slope increment and the rotation is approximated by a 
shearing transformation. The listing of CP (INF) from the equivalent free air 
solution gives a pressure distribution in the form of coefficients based on the 
corrected free stream conditions and is, of course, free of the influence of 
any nonuniformity in the wall-induced velocity perturbation. 

The final output gives the corrected lift coefficient and lists the 
result of transforming the experimental airfoil pressures to coefficients 
based on the corrected free stream conditions. The effect of any nonuniformity 
in the wall-induced velocity perturbations can be seen by comparing this 
distribution with the one resulting from the equivalent free air solution. 
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APPENDIX 


SAMPLE CASE 

The sample case consists of the application of program TWINTAN to a test 
case which was provided by M. Mokry and L. H. Ohman of the National Aeronautical 
Establishment, Ottawa, Ontario, Canada and has been used by them to compare 
several methods for Interference corrections. The experimental data were 
obtained in a test of the BGK-1 airfoil in the NAE blowdown wind tunnel with 
20 percent porosity top and bottom walls. The ratio of tunnel height to air- 
foil chord was 6. The airfoil was set at 2.56° incidence and tested at a 

£ 

reference Mach number of 0.784 and a Reynolds number of 21.03 x 10°. 

The input card images are listed and followed by the program output 
listing. The results of the model perturbation solution give a corrected 
Mach number of about 0.767 and show distributions of the wall-induced velocities 
VW and UW which are nearly uniform over the model chord. The results of the 
equivalent free air solution indicate an angle of attack correction of about 
-0.0112 radians which corresponds to -0.64°. If, however, the airfoil shape 
calculated in the tunnel flow solution is plotted and compared with the actual 
BGK-1 airfoil shape at 2.56° incidence, it is found that the correspondence 
between shapes is improved by adding downwash to the tunnel flow solution 
equivalent to a rotation of 0.25°. The corresponding angle of attack correction 
would then be -0.89°. 



INPUT FOR SAMPLE CASE 


SNPUT Jl*29, J3=51, J4=53, J5*75,M=34, DXMN*.0488, JU = 1, JD=75, KB*23 ,N*30, 
DYMN=.048333»KW*21, 

AMREF* .7e4, Cl* . 764 , C D «. 03, S l A -O ., SLR*0., JP-44, 

ICNSV-C, OPF(l) =2*1 .85,ITPC (1 ) -2,2,1$ 

$RDPX NW(l)=15,14,XW(l,l)*-7.85,-5.75,-3.35,-2.15,-1.25,-.65,.l,.7,1.15, 

1. 6, 2. 35, 2. 95, 3. 55, A. 15, A. 75, XW( 1,2) =-7.85,-5.75, -3. 3 5, -2. 15, -1.25, -.65, 
•1,.7,1.3,1.75,2.35,2.95,4.15>4.75,NM{1)«24,14,XM(1,1)*.025,.0402, .0602, 
.08015, .1002 ,.14,. 2, • 26, . 3402, . 400 4, . 4402 , . 5 004, . 56045, • 58035, .6002, .6201, 
.6501, .7001, . 750 3, .8003, . 8502 , . 9C0 4, . 95, 1 • » X M( 1, 2 ) * • 0198, • 05, . 1003, .1504, .2001 
.2501, .3496, .4 501,. 5 504, .65035, .7503, .8497, .95,1., 

CPW(1,1)»-. 0063, -.001 3, .0135, .0191, . 0064, -.004 9, - . 02 12 ,-. 0277, - • 0242 , - • 0245, 

-. 0 274, -. 0232, -. 0196, -. 01 72, -.0172, CP W( 1,2) *-.0064,. 01, .029, .04 22, .0 553, 

.0721, .0842 ,.0944, .0949, .0 90 5 » .0797, .0688 , . 039 6, . C257, C PM ( 1 , 1 ) 4455, 

-. 5 856, -.76e6,-. 9334, -.9919, -1.0309, -1.03 27, -1 .04 46, -1.0476, -1.0623 , -1 .082 8 , 
“1*112,— 1.1347,— 1.1416,— 1.1343,— .8942,— .5369,— .3707,— .2454,— .1454,— .044,. 0496, 

.1344, • 224, C PM (1,2) *.4226, . 1443 ,-. 0498, -. 1415, -. 1 906, -. 2088, -. 2049, -. 1172, 
-.0072, .1695, .3345, .4078, .3795, .224$ 



SAMPLE CASE OUTPUT 


SNPUT 

J1 

J3 

J4 

J 5 

H 

DXMN 

JU 

JO 

KB 

N 

DYMN 

KW 

AMREF 

CL 

CO 

SLA 

SLB 

JP 

ICNSV 

RFM 

RFA 

ORF 

OXT 

ITRC 

EPS 

I MX 


* 29# 

- 51# 

* 53# 

- 75# 

« 34# 

« • 488E-01# 

- 1 # 

« 75# 

- 23# 

* 30# 

- *4 8333E-01# 

- 21 # 

* .784E+00# 

» .764E+0C# 

- .3E-01# 

■ 0 . 0 # 

= 0 . 0 # 

* 44# 

- 0 # 

» • 2 E +00# 

= • 5 E + CO# 

« .185E+01# • 1 85 E +01 # .18E + C1# 
■ • 5E+00# .5E+00# • 5 E +00# 

- 2 > 2 » 1 # 

* • 5E-0 5# .5E-05# .5E-05# 

- 350# 500# 200# 


SEND 



J 

X(J) 

FKJ) 

F2{ J) 

XMID(J) 

1 

-• 78 1 65 E +01 

.63815E+00 

. 74894 E+00 

-.71489E+01 

2 

-.64734F+01 

. 86860E + 00 

•99711E+00 

59719E+01 

3 

— • 5 4660 E + 01 

•11345E+01 

•12807E+01 

-.50756E+01 

4 

-.46825E+01 

• 14358E+01 

•15998E+01 

-. 43700E+01 

5 

-.40557E+01 

• 1 7726E+01 

• 19543E+01 

— • 37998E+01 

6 

-•35426E+01 

•21449E+01 

•23443E+01 

-• 33296E +01 

7 

-.31155E+01 

•25526E+01 

• 27698E+01 

-.29349E+01 

8 

-• 275 38 E +01 

• 29958E+01 

.32306E+01 

-. 25991E+01 

9 

-.24439E+01 

•34744E+01 

• 37270E+01 

-• 23097E+01 

10 

-.21753E+01 

. 39885E+01 

•42588E+01 

-.20578E+01 

11 

-.19402E+01 

• 453 60E + 01 

• 48 260E+01 

-.18366E+01 

12 

-.17328E+01 

• 51230E+01 

.542 87E + 01 

-.16407E+01 

13 

-.15484E+01 

•57434E+01 

•60669E+01 

— • 14660E + 01 

14 

-.13835E+01 

•63993E+01 

•67405E+01 

-. 13093E+01 

15 

-.1235CE+01 

•70906E+01 

• 74496E+01 

-.11679E+01 

16 

-•1 1007 E+01 

. 78174E+01 

, 81 941E + 0 1 

10397E +01 

17 

-.97862E+00 

•85796E+01 

* 89740E+01 

-.92290E+00 

18 

-• 8671 3 E +00 

•93773E+01 

• 97 894E + 01 

-• 81606E +00 

19 

-•76493E+00 

•10210E+02 

• 10640E+02 

-•71794E+00 

20 

-.67091 E+00 

. 11079E+02 

•11527E+02 

-•62753E +00 

21 

-•58412E+00 

•11963E+02 

•12448E+02 . 

-• 54396E+00 

22 

-.50376E+00 

♦12923F+02 

•13406E+02 

46647E +00 

23 

-.42914E+00 

• 13 8 98 E +02 

• 14398 E + 02 

-.39442E+00 

24 

359 67 E +00 

. 1 4908 E + 02 

•15426E+02 

-. 32726E+00 

25 

-.29483E+C0 

•1 5954E+02 

« 16490E+02 

-•26451E+00 

26 

-• 23417 E +00 

•17G35E+02 

*175 89E + 02 

-. 20574 E+00 

27 

-• 1773GE + 00 

•16152E+02 

• 18724E+02 

- . 1 5060E +00 

28 

-• 12388 E+00 

. 19 3C4E + 02 

.1 9894E+02 

-.98745E-01 

29 

-.7360GE-C1 

. 20492E+02 

.20492E+C2 

- *4 9200E— 0 1 

30 

— • 24 800 E— 01 

.204926+02 

.204926+02 

-.400006-03 

31 

•24000E-01 

•20492E+02 

• 20492 E+02 

. 4 8400E-01 

32 

.728 00 E-0 1 

.20492E+02 

•20492E+02 

. 97200E-01 

33 

.1216CE+00 

•2C492E+02 

* 20492E+02 

. 14600E+00 

34 

. 1 7040 E +00 

. 2C492E+02 

.20492E+02 

. 19480E+00 

35 

• 2 192C E + 00 

.20492E+02 

• 20492 E+02 

• 24360E + 00 

36 

« 26800 E +00 

. 20 49 2 E+02 

•20492E+02 

• 29240E+00 

37 

• 3 16 60 E +00 

•20492E+02 

• 20492E+02 

* 34120E + 00 

38 

•3656GE+00 

.2C492E+02 

•20492E+02 

• 39000E +00 

39 

• 41 4 40 E +00 

•20492E+02 

* 20492E+02 

. 43880E+00 

40 

.4 63 20 E+00 

♦20492E+02 

•20492E+02 

• 48760E+00 

41 

. 5 12C0E+00 

. 20 4 92 E +02 

. 20492E+02 

• 5 3640E + 00 

42 

. 560 80 E +00 

•20492E+02 

,2049?E+0? 

• 58520E+00 

43 

•6096CE+00 

• 2C 4 92 E+02 

• 20492E+02 

• 6 3400E + 00 

44 

. 65840E+00 

• 20492E+02 

.20492E+02 

. 68280E+00 

45 

• 7072C E + 00 

•20492E+02 

• 20492E+02 

. 73160E + 00 

46 

.75600 E +00 

.20 49 2 E+02 

.20492E+02 

. 78040E +00 

47 

. G048CE+00 

« 204 92 E + 02 

•20492E+02 

. 82920E +00 

48 

.8 53 60 E+00 

• 20 49 2E + 02 

•20492E+02 

. 87800F+00 

49 

•9024CE+00 

•2C492E+02 

* 20492E+02 

,92680.E + 00 

50 

•95120E+00 

.2C492E+02 

»20492E+02 

• 97560E +00 
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51 

• 1O0C0E+01 

•20492E+02 

• 20492 E+02 

•10244E+01 

52 

• 104 88 E +01 

• 20492E+02 

•20492E+02 

. 10732 E +01 

53 

• 1 0976E+01 

•2C492E+02 

* 19394E+02 

•11227E+01 

54 

•11479E+01 

•193C4E+02 

, 1 8724E + 02 

. 11746E+01 

55 

• 12013E+01 

•18152E+02 

•17589E+02 

• 12297E+01 

56 

• 12582E+01 

•17035E+02 

, 1 64 90E + 02 

• 1 2 885 E +01 

57 

• 131 88E+01 

•159 5 4E +02 

• 15426E+02 

• 135 12E + 01 

58 

• 1 38 37 E+01' 

•14908E+02 

•14398E+02 

• 1 41 84 E+01 

59 

• 145 31 E + 01 

• 13 898E+02 

•1 3406E+02 

• 14904E+01 

60 

•15278E+01 

•12923E+02 

• 1 244 8E + 02 

. 15679E+01 

61 

• 160 81E + C1 

•11983E+02 

• 11527E+ 02 

• 16515E +0 1 

62 

•1 6949E+01 

• 11079 E + 02 

•10640E+02 

• 17419E+01 

63 

• 178 69 E +01 

•10210E+02 

• 97894 E+01 

• 18400E+01 

64 

•18911 E+01 

•93773E+01 

• 89740E+01 

• 19468E+01 

65 

• 20C26E+01 

• 857 96E+01 

• 81941E+01 

• 20636E +01 

66 

• 21 247 E+01 

•78174E+01 

,744 96E + 01 

. 21918 E+01 

67 

•22590E+01 

• 70906E+01 

•67405E+01 

• 23332E+01 

68 

• 24075E+01 

•63993E+01 

• 60669E + 01 • 

•24899E+01 

69 

• 25724E+01 

•57434E+01 

• 54287 E+01 

• 26645E+01 

70 

• 27 568 E +01 

. 51230E+C1 

•48260E+01 

• 28604 E +01 

71 

• 2 9642 E +01 

•45380E+01 

• 42588 E + 01 

• 308 16E + 01 

72 

• 31 993 E +01 

. 3988 5 E +01 

•37270E+01 

• 33334E+01 

73 

• 34679E+01 

•34744E+01 

• 32306E+01 

•36227E+01 

74 

• 37778 E+01 

•29958E+01 

•27698E+01 

. 39584E+01 

75 

•41395E+01 

• 25 526 E + 01 

• 23443 E +01 

• 43527E+01 


K 

Y( K) 

G1(K) 

G2(K) 

1 

0. 

•2C690E+02 

•20006E+02 

2 

•50C00E-01 

•19333E+02 

• 1 8 67 3 E+02 

3 

• 103 57 E +00 

•18023E+02 

• 173 85E + 02 

4 

•16111E+0C 

•16759E+02 

. 161 44E + 02 

5 

. 22308 E + 00 

. 155 40E+02 

• 14948E + 02 

6 

• 29000E + GC 

•14368E+02 

• 1 3799 E + 02 

7 

. 362 50 E +00 

•13241E+02 

•12695E+02 

8 

•4413CE+00 

•12 161E+02 

• ,11638 E + 02 

9 

. 52727E+00 

•11 127E+02 

•10627E+02 

10 

• 62142 E + OC 

•10138E+C2 

• 966 10E + 01 

11 

• 72 500 E+00 

•91955E+01 

• 87414E +01 

12 

• 83947E + CC 

• 62989E+01 

• 78679E + 01 

13 

•96666c +00 

• 74483E+01 

•70403E+01 

14 

• 11C88E+01 

•66437E+01 

• 6 2 5 87 E +01 

15 

. 1 2 6 87 E +0 1 

•5P851E+01 

•55230E+01 

16 

•1450CE+01 

•51724E+01 

• 48 3 34 E + 01 

17 

•1 6571 E+01 

• 45058E+01 

•41897E+01 

18 

• 1 8961 E + 01 

•38651 E+01 

• 3 59 20 E +01 

19 

• 217 50E + 01 

. 23104E+01 

•30403E+01 

20 

• 2 5045 E+01 

•27816E+01 

• 25 345E+01 

21 

• 29000 E+01 

• 22989E+01 

•20747E+01 

22 

•33833E+01 

•16621E+01 

• 1 6609 E+01 

23 

•39875E+01 

• 14713E+01 

•12931E+01 
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AIPFOIl CP INTERPOLATED AT MESH MIDPOINTS 


-.66018 

-1.13080 

-.98743 

-1.15017 

-1.03323 

-.70513 

-1.03281 

-.39726 

-1.04092 

-.29534 

-1.04717 

-.18227 

-1.04767 

-.08718 

-1.05809 

.00947 

-1.08208 

.09477 

-1.10536 

.17975 

.15611 

-♦02606 

-.04324 

.04798 

-.13498 

.13809 

-.18724 

.23050 

-.20728 

.31041 

-.21550 

.36387 

-.20857 

•39603 

-.17643 

.42325 

-.12863 

.41531 

-.08033 

.30876 


I 

JMX 

KMX 

ISUP 

NS 

DPHIMX 

TTE 

10 

5 8 

-1 

0 

83 

* 14972 E + 00 

0. 

20 

59 

-1 

0 

98 

•32637E-01 

0. 

30 

60 

-1 

0 

104 

•5791 8E-02 

0. 

40 

43 

10 

0 

.08 

•2822 5E -02 

0. 

50 

27 

-12 

0 

94 

-.10169E-02 

0. 

60 

31 

3 

0 

97 

-•23 61 9E -02 

— .1011 IE +00 

70 

31 

3 

0 

96 

-• 10462E-02 

— •10111E + 00 

BO 

31 

3 

0 

95 

- .9900 1 E-03 

-.40739E-01 

90 

29 

1 

0 

95 

— • 10775 E— 02 

— *36273 E-03 

ICO 

29 

1 

0 

94 

-• 691 97 E-03 

• 1 6 329 £-01 

110 

27 

-9 

0 

96 

-.35498E-03 

* 14807E— 01 

120 

28 

-10 

0 

96 

-•1 9063E-03 • 

1 5835E-01 

130 

25 

-13 

0 

96 

-.94511 E-04 

•1 5230E-01 

140 

2 5 

-11 

D 

96 

-•47021 E— 04 

. • 15148E — 01. 

150 

66 

-3 

0 

96 

— • 13578E-03 

. 1 5 0 8 8 E —0 1 

160 

57 

8 

0 

96 

-•4 273 9 E-04 

. 1 4991E— 01 

170 

47 

-10 

0 

96 

— • 23 393 E-04 

•14965E— 01 

1B0 

38 

-14 

0 

96 

-.12391E-04 

•15014E— 01 

190 

2 9 

-15 

0 

96 

— .62766E-05 

•1 5024E-01 

195 

24 

-13 

0 

96 

- *481 87 E— 05 

•15020E-01 



PESUITS IN TUNNCL FLOW 


CPSTAR « <> 81 83E *00 


X( J-.5) 

YMOD 

CP(REF) 

U 

J 

X(J) 

V 

DELTA PHI 

DELTA V 

-.00040 

0.00000 

.81668 

-.44457 

31 

. 024C0 

•64698 




o.oocoo 

1.04505 

-.69171 



-.61341 

.01206 

1.26040 

.04640 

.03157 

-.6553 1 

.43819 

32 

.07280 

.23048 




-.02993 

.14630 

.00250 



-.21816 

.03332 

.44863 

.09720 

.04282 

-.98835 

.50106 

33 

.12160 

.10789 




-.04050 

-.04325 

.03945 



-.15799 

.05505 

.26588 

.14600 

• 04 RC 6 

-1.03534 

.50885 

34 

.17040 

.05745 




-.04829 

-.13500 

.0764? 



-.11232 

.07695 

.16978 

.19480 

.05089 

-1.03288 

.50379 

35 

.21920 

.02607 




-.05377 

-.18725 

.09705 



-.07635 

.09680 

.10242 

.24360 

.05216 

-1.04089 

.50665 

36 

.26800 

-.00159 




-.05750 

-.20729 

.10406 



-.04819 

.11645 

.04660 

.29240 

.05208 

-1 .04715 

.50982 

37 

.31680 

-.02213 




-.05985 

-.21550 

.10660 



-.02228 

.13612 

.00015 

. 34120 

.05100 

-1.04768 

.51066 

38 

.36560 

-.03753 




-.06094 

-.20857 

.10260 



.00332 

.1560^4 

-.04085 

.39000 

.04917 

-1.05821 

.51697 

39 

.41440 

-.05530 




-.06077 

-.17643 

.00705 



.02563 

.17702 

-.08093 

.43880 

.04647 

-1.C8205 

.53085 • 

40 

•46320 

-.07704 




-• C5952 

-.12863 

.06424 



.04327 

.19979 

-.12031 

.48760 

.04271 

-1.10508 

.54521 

41 

.51200 

-.09771 




-.05741 

-.08033 

.04126 



.05997 

.22438 

-.15769 

.53640 

.03795 

-1.13017 

.56146 

42 

•56080 

-.11982 




-.05449 

-.02606 

.01540 



.07776 

.25103 

-.19758 

. 58520 

. C32 1 C 

-1.14320 

.57180 

43 

.60960 

-.13630 




-.05069 

.04798 

- .02060 



.09196 

.27994 

-.22827 

.63400 

.02545 

-.70622 

.34367 

44 

.65840 

.02360 




-.04620 

.13808 ' 

-.06594 



.09698 

.29993 

-.07338 

.68280 

.02660 

-.40090 

.19637 

45 

.707 20 

-.11440 




-.04147 

, .23046 

-.11457 



.08910 

.31510 

-.20350 

.73160 

.02102 

-.29546 

.15453 

46 

.75600 

-.15059 




— • C 3 71 2 

.31037 

-.15875 



.06730 

.33039 

-.21790 

.78040 

.01367 

-.18235 

.10348 

47 

.80480 • 

-.16715 




— • 03 3R4 

.36385 

-.1*8963 



.03916 

.34469 

-.20631 

.82920 

.00551 

-.08720 

.05860 

48 

.85360 

-.17707 




-.03193 

.39602 

-.20857 



.01038 

.35773 

-.18744 

.87800 

-.00313 

.00947 

.01076 

49 

.90240 

-.17726 




-.03142 

.42325 

-.22424 



-.03276 

.36920 

-.14450 

.92680 

-.01178 

.09476 

-.03350 

50 

.95120 

-.17127 




-.03302 

.41527 

-.21831 



-.08807 

.37822 

-.08320 

.97560 

-.02014 

.17971 

-.07975 

51 

1.00000 

-.15391 




-.03732 

.30869 

-.15728 



-.08451 

. .36200 

-.06940 

1.02440 

-.02765 

. 








ro -.04144 



OUTER BOUNDARY VELOCITIES 



X - 

- *71 4 89E +01 

YU 

YV 

U 

0. 

• 24993F-01 

-•14172E-02 

•50000E-01 

•76777E-01 

-*1 3 974E-02 

* 10 3 57E + 00 

• 1 3233 E +00 

-.13740E-02 

• 1 61 1 1 E + 00 

* 1 9208E +00 

-.13463E-02 

• 223 OGC + OO 

• 2 5652 F +00 

-•131 36E-02 

• 2 9000E + 00 

.32623E+00 

-.12749E-02 

• 3 62 50 F + 00 

• A 01 88 E +00 

— •12 29 IE- 02 

• 4413CE+00 

• 4 8426E +00 

-.11745E-02 

•52727E+C0 

•57432E+00 

-•11096E-02 

• 62142E+00 

• 6731 8F +00 

— • 1 0320E-02 

•7250CE+00 

•78219E+C0 

-.93881E-03 

• 83947E + G0 

• 90302 E + 00 

— .8 2653 E-03 

•96666F+00 

• 10377E + 01 

-•69052E-03 

•11088E+01 

•11P87E+C1 

-.52485E-03 

•12687E+01 

•13593E+C1 

— • 3 21 81 F-03 

•14500E+01 

.1 5534E+01 

-•713 63 E -04 

• 16571 E + Gl 

• 1 7765 E +01 

• 2 3 9 66 E-0 3 

* 1 8961 E + C1 

• 2 0? 5 3F +01 

•62847E-03 

• 2 1750E + C1 

• 2 3 3 94 F +01 

• 1 1 1 7 3 E-02 

•25045E+01 

• 270 18E + 01 

*17 3 3 5 E-02 

0. 

-.24993E-01 

-•1 4172 E-02 

- • 50000 E-01 

-.76777F-01 

-.14351 F-02 

-.10357E+00 

— •1 3233E +C0 

-• 14 520E-02 

-• 16111E+00 

-.1920FE+00 

-•1 4 677 E-02 

-* 22308 E + 00 

-.25652F+00 

- .1 4 81 6 E-02 

-• 290Q0E + CO 

-.32623E+00 

-•14933 F-02 

-* 36250E+00 

— .401888 + 00 

-*1 5020 E-02 

-•44130E+00 

— • 4 0426 E +00 

— • 1 5069 E-02 

-• 52727E+00 

-.57432F+C0 

-*1 5065 E-02 

-•62142E+00 

— •6731GE+00 

-•14997E-02 

-.72500E+00 

-.78219E+00 

-•14 8 44 E-0 2 

-•8 39 47E+00 

-•902C2r+00 

-.145 82 E-02 

-• 96666E + 00 

-.10377E+01 

— •14 1 77 E-02 

-♦ 1 108 8 E +01 

-.11887F+G1 

-*13 5 86 E-02 

-•12687E+01 

— •13593E+01 

-• 12753E-02 

-•145C0F+01 

-.15534F+C1 

- • 1 1 602 E-0 2 

-.16571E+01 

-.17765E+01 

-•10033 E-02 

-• 1 89 61E + 0 1 

'-.20353E+01 

-.79135E-03 

-• 21 750 E +01 

- • 2 3 394E +01 

- • 506 90 E-03 

-•25045E+01 

-•27018E+01 

— •12016E-O3 


— • 70 16 5 E + 01 

• 395 84 E +0 1 

•41395E+01 

V 

U 

V 

• 63 647 E-03 

-.71975E-02 

-.14018E-01 

• 61 3 60 E-03 

-.69205E-02 

-.14145E-01 

• 58 868 E-03 

-• 66230E-02 

— • 14277E-01 

• 56 15GE-03 

-.63027E-02 

— .14A15E-01 

•531 80E-03 

— • 59569E— 02 

-.14559E-01 

. 499 32 E-03 

-.55827E-02 

-.14710E-01 

• 46 380 E-03 

-.51764E-02 

— «14867E-01 

• 424 92 E-03 

— • 47340 E-02 

— • 15030E-01 

• 3 8 24 0 E-03 

-.42507E-02 

-.15200E-01 

• 33 5 95 E-03 

-•37208E-02 

— • 15 3 75E-01 

• 2 8 5 3 7 E-03 

-•31378E-02 

- • 1 5 555E-01 

• 23056 E-03 

— • 249 37 E-02 

-.15739E-01 

• 1 71 67 E-03 

— .1 7791 E-02 

-.15923E-01 

• 10 92 8E-03 

-•98259E-03 

-.16105E-01 

• 447 96E-04 

— • 90493 E-04 

-.16278E-01 

-* 1 9009 E-04 

• 914 2 2 E — 03 

-.1 6434 E-01 

-.76P23E-04 

•20529E— 02 

-.16559E-01 

-• 118 7 6 E-03 

• 33530E-02 

-•16632E-01 

— •12 641 E — 03 . 

• 4 8 523 E-02 

-.16616E-01 

— • 65014 E-04 

• 66061E-02 

— • 1 64 56E-01 

• 6 5 81 0 E-03 

-•71975E-02 

- »13 894E-01 

• 6 6026 E-03 

-•74743E-02 

-.13762E-01 

• 70341 E-03 

— *77700 E-02 

-.13617E-01 

• 72764E-03 

— • 808 66 E— 02 

— • 13456E— 01 

• 75 29 2E-03 

- • 84262 E-02 

-.13279E-01 

• 7791 6 E— 03 

-.67913E-02 

-.13082E-01 

• 80624 E-03 

— • 91 850 E— 02 

— • 12861E— 01 

• 83391 E-03 

— • 96105E-02 

-.12614E-01 

♦ 86 179 E-03 

-.10072E-01 

-.12334E-01 

• 88 92 9E-03 

— • 105 74E-01 

-.12016E-01 

• 91 5 4 8E-03 

-•11 121E-01 

-• 11652E-01 

• 93897 E-03 

-•11722E-01 

-.11232E-01 

• 95765 E-03 

— •12382 E-01 . 

-.10744E-01 

• 968 27 E-03 

- *13113 E-01 

-• 10171 E-01 

• 96582E-03 

- .13 92 8 E-01 

-.94920E-02 

• 9424 8 E— 03 

-.14842E-01 

— • 86774E— 02 

* 09 5 78 E-03 

- • 15 880E-01 

-.76869E-02 

• 77547 E-03 

— • 17074 E— 01 

-.64636E-02 

• 57767 E-03 

- • 1 8475E— 01 

-.49254E-02 

.•233 88E-03 

— • 20168 E-01 

-.29493E-02 




Y « 

•29000E+01 

•29000E+01 

-.29000E+01 

-•29000E+01 

XU 

XV 

U 

V 

U 

V 

71499E + 01 

-•78165E+01 

• 2 5072E-02 


•36941E-03 


59719E+01 

— • 64734E +01 

• 10 247E-02 

• 265 50E-02 

-.41939E-02 

— .17655E— 02 

50756E+01 

— #54660 E +01 

“.80003E-03 

• 46979E— 02 

-.73929E-02 

-.16798E-02 

43700 E +01 

-.46825E+C1 

“.28168E-02 

•64776E-02 

— .99856E-02 

— .1494 8E-02 

3799 PE + 01 

-•40557E+01 

“• 48768E—02 

• 7991 2 E-02 

-.12365E-01 

— .14866E— 02 

33296E+01 

-•35426E+01 

-.68958E-02 

•92090E— 02 

-.14649E-01 

— • 16209E— 02 

29349E+01 

-•31 1 55 E +01 

“• 867 39F-02 

• 100 31 E-01 

— •16818E-01 

-.18286E-02 

25991E+01 

-.27538E+01 

-• 973 69E-02 

• 10334 E— 01 

— #1 8734 E-01 

— » 2050 3E— 02 

23097E+C1 

— *24 439E + C1 

“•99262 E-02 

• 10 31 3 E-01 

— .20346E-01 

— .23374E— 02 

2057PE+01 

— #21753E +01 

“• 9 2929 E“0 2 

• 1 01 89E— 0 1 

-.21663E-01 

— .27598E— 02 

18 36fcE +01 

-•19402E+01 

-.80934E-02 

• 1 01 79 E— 01 

-.22837E-01 

-.34230E-02 

1 6 4 07E + 01 

-•1 7328E+01 

“• 66636E-02 

• 10391 E— 01 

-.24030E-01 

— .43724E— 02 

1 4660 E +01 

-.15484E+01 

“• 52005 E“02 

• 10758 E-01 

-.25495E-01 

— • 55 152E— 02 

13093E+01 

- *1 3 P3 5 E +01 

“.38100E-02 

• 11 19 3 E— 01 

-.27112E-01 

-.67484E-02 

11679E+01 

-•12350E+01 

“• 25 4 27E-02 

• 11614 E — 01 

-.28919E-01 

-. 79746E-02 

10397E+C1 

-.11CC7F+01 

-•13881 E-02 

• 11936 E-01 

— • 3078 7E-01 

-•90428E— 02 

92290E+00 

- • 978 62 E +00 

-.31399E-03 

• 12 1 1 3E— 01 

-.32560E-01 

— .98811E— 02 

8 1 60 f E + 00 

-•06713E+CO 

•701 50E-03 

•12131 E-01 

— • 34 15 IE-01 

—.1051 9E— 01 

71 7 94 E +00 

— #7649 3E +00 

• 16726E-02 

• 1 1 994 E— 01 

— • 3551 5 E-01 

— .10996E— 01 

62753E + 00 

-.67091E +00 

• 2 6084 E-02 

. 1 1 70 7E — 01 

— *3663 6 E-01 

— . 11353E— 01 

54396E+00 

-.58412E+00 

• 3 5 100E-02 

• 11 2 8 8 E-01 

— .37534 E-01 

-.1164 9E— 0 1 

46647 E +00 

-.50376E+00 

•43696E-02 

• 1075 7 E — 01 

— .38267E-01 

-• 11951E— 01 

39442E+00 

-.42914E+00 

• 5 1 82 1 E-02 

• 101 27E-01 

— .38879E-01 

-.12268E— 01 

32726E+00 

-*3 5967E + C0 

• 59454E-02 

• 9406 1 E— 02 

— .39404E-01 

— .1 2599E-01 

26451E+00 

-.29483E+00 

• 66 5 88 E-02 

• 8606 9E-02 

— .39867E-01 

— .12939E-01 

20574E+00 

- • 2 3417 E +00 

• 73 233E-02 

• 77402 E-02 

— .40287E-01 

— .13284E-01 

1 50 60 E +00 

-•1773OE+C0 

• 79403 E-02 

• 68 173E-02 

— .40675E-01 

-•13630E— 01 

96745E-01 

-#1238eE+00 

• 8 5 1 18 E-02 

• 5 84 9 2 E-02 

— .41043E-01 

-.13971 E— 01 

49200E-01 

— •7360CE-C1 

• 90 477 E-02 

• 4 9348E— 02 

-.41402E-01 

-.14316E-01 

40000E-03 

- #24 PCGE-01 

• 9563 5E-02 

• 379 1 5E-02 

— .41768E-01 

— .14641E— 01 

4 8 4 OOE-Ol 

• 24000 E-0 1 

• 10065E-01 

• 26525 E-02 

— .42151 E-01 

-• 14984E— 01 

97200E-01 

• 72 8 OCE -01 

• 10 550 E-01 

• 14334 E-02 

— .42557E-01 

-•15332E-01 

14600E+00 

• 1 2 1 60 E + 00 

• 1 1017E— 01 

• 13 083E-O3 

— • 4298 9E-01 

-•15678E-01 

1 94 80 E +00 

• 1704CE + 00 

• 1 1 460E-01 

-• 12 2 30E-02 

— .43442E-01 

-.16005E-01 

2436CE+00 

• 21 920 E +00 

• 11 8 76E-01 

— • 2643 6E-02 

— . 43910E-01 

-.16304E-01 

292^0E+00 

.26eC0E+00 

• 1 2 260E-01 

— •4114 5 E-02 

— • 44386E-01 

-. 16574E— 01 



•34120E+C0 
•39000E+00 
•43880E+00 
•48760E+C0 
• 53640E+00 

• 5 8 5 20 E + 00 

• 6 3400 E + 00 
•68280E+C0 

• 73 160E + 00 

• 78 040 E + 00 

• 82 92 0E + 00 

• 8 7800 E + 00 

• 926 80 E + 00 

• 9 75 6 OE + 00 
•10244E+G1 

• 1073 2 E + 01 
•11227E+01 

• 1 1746E + 01 
•12297E+01 
•12885 E +01 

• 13 51 2 E + 01 

• 141 84E + 01 
•14904E+C1 
.15679E+C1 
•16515E+C1 
•17419E+01 

• 1 84 00 E + 0 1 
•I9468F+01 

• 2063 6 E +01 
•21910E+C1 
•23332E+01 
•24899E+01 

• 266 A 5 E +01 

• 2 660A E + 01 

• 300 16 E +01 

• 3333 A E +01 
• 36227E+01 

• 3958A E +01 


•3168CF+00 

•3656CF+00 

•A1AACF+C0 

• A 63 20 E +00 

• 5 1 200F + 00 

• 560 60 F +00 
•60960F+00 
•656ACF+00 

• 7072CE +00 

• 75 6 0CF + 00 
•80A8CE+00 
•8536CF+00 

• 902 A 0 E +00 
•9512CE+00 

• 1 OCCCE +C1 
•10A6GE+C1 
« 10976F+01 

• 1 1 A 79E +0 1 
•12013E+01 

• 1 2 5 02 E +01 
•13188E+C1 

• 13837E + 01 

• 1A531 E +0 1 
•15276E+01 

• 1 6081 E + 01 
•16949E+C1 
•17809E+O1 
.1 891 1 E +01 
•20026E+C1 

• ?1 2A7E +0 1 
•2259CE+01 
•2A075E+01 

• 2572AE + 01 
•27566E+01 
•29642E+01 

• 3199.2F + C1 
•34679E+01 

• 3 7778E + 01 


•12608E-01 
• 1 2916E-01 

• 1 3179E-01 
. 13393E-01 

• 1 3 5 55 E-01 
. 13660E-01 

• 1 3 705 E-01 

• 1 3 6 86 F-01 

• 1 3 599 F-0 1 
•13451E-01 
• 1 3255 E-01 
•13024E-01 

• 12 770E-01 

• 1 2 507E-01 

• 122 A7F-01 

• 12 003E-01 

• 1 1 7 8 5 E— 0 1 
•11603E-01 

• 1 1 468 E-01 
•11384E-01 

• 11 356 E-01 

• 1 13 87 E-01 

• 1 1 A78 E-01 

• 1 1 631 E-01 
•11840E-01 

• 12099E-01' 

• 1 2 396 E-01 
• 1 2 710 E-01 
•1300AE-01 

• 13219E-01 
•1 3260E-01 
•13004E-01 
•12A36E-01 
•* 116AAE-01 
•10832E-01 
*1011 5E-01 
.9A160E-02 
• 87079E-02 


— • 562A 6E-02 
-•71620E-02 

— • 8 71 A0E-02 
-•10267E-01 

— •11 807 E-01 
-•13322E-01 

— • 1 A798 E-01 
— • 16226E-01 

— • 1 75 99E-01 

18923E-01 

— • 20205 E-01 
— • 21 AA0E-01 

— •22 627 E-01 
-.23762E-01 

— • 2 A8 A A E-01 

— •25 G71E-01 

— • 268 39 E-01 

— •27789E— 01 

— • 2 871 1 E-01 

— •2 95 82 E-01 

— • 30 3 85 E— 01 
•-• 31100E-01 

-•31709E-01 
— • 32188E— 01 
-.3251AE-01 
—•32 66 A E-01 

— • 3261 3 E-01 

— • 32333 E— 01 
— • 31 796E-01 

— • 30 977E-0 1 

— • 2985AE-01 
-.28A38E-01 

— • 268 A 5 E-01 

— •251 85 E-01 
-• 235A2 E-01 

— • 218 73 E-01 

— ♦200 85 E— 01 
• 18 217E— 01 


— • AA865 E-01 

— • A5338E— 01 
-•A5801E-01 

— • A 62 A6E— 01 

— • A6667 E— 01 

— • A705 7E-01 
— • 47411 E-01 

— • A7720E-01 
-•A79Q0E-O1 
- • A 81 90 E — 01 

— • A 83 52E-01 
- .A 8A7 0 E-01 
-• A85A7E-01 

— • A8 58 A E— 01 
— • A8586E— 01 

— • A 8 5 56E-01 

— • A 8A 95E — 01 

— • A 8A01 E— 01 

— • A 8272E-01 

— •A8106E— 01 
-• A7902E— 01 
-.A7651E-01 
-.A73A1E-01 

— • A6955 E— 01 
-.A6A73E-01 

— • A 5867 £-01 
- • A5 105 E— 01 
-.AA173E-01 

— • A3071 E-01 
- • A 1 81 2 E— 01 
- • AOA32E-01 
— • 38979E— 01 

— • 37A09E — 01 
-•3 55 87E-01 

— • 3 3289 E— 0 1 
-•3033AE-01 

— * 26665 E— 01 

— •22303E-01 


-•168 11E-01 

- • 170 1 AE — 01 
-•17183E-01 
-• 17316E-01 
-• 17A1 3E-01 

- • 17 A75 E— 0 1 
-•1750 AE — 01 
-•17A99E-01 

- • 17A6A E— 01 
- • 1 7 A 1 2 E — 0 1 
-• 1735 2E-01 
-• 1 7285E-01 

- • 172 1 AE-0 1 
- • 1 71 39E— 01 
-• 1 7062 E-01 
-• 1698 AE— 01 
—•1 6901 E-0 1 
-• 16817E— 01 
-• 16725E— 01 
-•16623E— 01 
-• 16509E— 01 
-• 1 6368E— 01 
— •161 89E— 01 
-.1$965E-01 
-• 156 86E-01 

- • 153 A2E-01 

- • 1A930E-01 
-• 14A71E-01 
-•13975E-01 
- • 1 3 A A5 E— 01 
-• 12 8 0 IE-01 
-• 12258E— 01 
- • 11 A63E— 01 

- • 10A01E-01 

- • 9015 4E-0 2 
-•73698E— 02 
-•55885 E— 02 
-• 37565E— 02 



PHI IN VICINITY OF AIRFOIL 


J 

PHI ( J#3>2) 

PHI ( J # 2# 2 ) 

PH I ( J * 1 * 2 ) 

PHI ( J, 1,1) 

PHI( J#2#l) 

PHI(J,3#1) 

29 

-•14693 F + 00 

-•14677E+00 

-•14428E+00 

- • 1442 0 E + 00 

— • 13808E + 00 

-• 13021 E+00 

30 

-• 15732E+00 

-•1 601 0 E +00 

— • 15975E +00 

-•1 597 5 E+00 

-. 14830E+00 

-.13661E+00 

31 

- • 1 679 5 E +00 

— • 17693E +00 

-•193 5 IE +00 

— • 18145E + 00 

-•15778E+00 

— • 14087E + 00 

32 

-•17422E+00 

-•ie281E+00 

- • 1 933 9E +00 

- • 1 6007E +00 

-• 148 55 E +00 

-• 13620E + 00 

33 

— • 1768 3E + 00 

-•18384E+00 

— • 191 46E +00 

-•13561E+00 

-.12948E+00 

-•12171 E+00 

34 

-•17690F+00 

— • 1 8 2 2 5E +00 

-•18773E+00 

-•11078 E + 00 

- • 10780E+00 

-•1035 2E +00 

35 

-• 17539E+00 

-•17922F+00 

— • 18300E+00 

— • 86198 E — 01 

— • 84965E-01 

-• 83047E-01 

36 

-•17303E+00 

— • 1 75 5 2E +00 

— • 17792E +00 

— • 61473E-01 

-• 615 13 E-0 1 

— •61214E-01 

37 

-•17032E+00 

— •1715 8E + 00 

-•17272E +00 

— •36 594 E-01 

-•37656E-01 

— • 3 85 56E-01 

38 

- • 16767E +00 

-•16778F+00 

-•1 6771 E +00 

- • 1 1 674E-01 

— •13538E-01 

— •15 379E-01 

39 

-.16547F+Q0 

-•16459 F +00 

-• 2 6346F +00 

• 1 3 5 5 4 E-0 1 

• 10880E-01 

• 820 1 IE-02 

40 

-•164C4F+00 

-•16231F+00 

-• 1 6033 E +00 

• 3 9 46 0E-01 

• 3 5 81 5 E-01 

•32240E-01 

41 

-.16353F+00 

— •16 1 09 E + 00 

-•15 831 E +00 

« 66066 E-01 

• 61406E-01 

• 56 854 E-01 

42 

- • 164 1 5 F +00 

— •1 61 1 OE +00 

-•15756E+00 

• 9 3 46 5 E-01 

• 87744 E-01 

• 82 144E-01 

43 

-.16607E+00 

-•16265E +00 

-.15857E+00 

•12137E+00 

•11473E+0O 

• 10811E+00 

44 

-.16944E+00 

- *1 6601 E +00 

— • 1 6179 E+00 

• 1 38 1 4 E+00 

. 13677E+00 

• 132 17E + 00 

45 

-•1742 3 F + 00 

-•17119E+00 

— • 1 6738 E +00 

• 1 4772 E +00 

•14211E+00 

•13598E+00 

46 

-•18026E+00 

-• 17801 F +00 

- • 1751 2 E + 00 

• 1 5 526 E + 00 

• 14797E +00 

• 14057E +00 

47 

— • 18 72 4 E + 00 

- *1 860 3 E +00 

-•18438E +00 

•16031E+00 

•15225E+00 

• 14421 E+00 

48 

-•19481 E + 00 

—•1948 1 E+00 

— •1 9456E +00 

• 16317 E+00 

• 15472E +00 

• 1 4642E +00 

49 

— • 20260E+00 

-.20396E+00 

— • 20550E +00 

• 16370E + CO 

•15529E+00 

• 147 10E + 00 

50 

-• 21 010E+00 

— • 21 273E+00 

-• 2161 5E +00 

• 16206 E + 00 

• 15405 E +00 

• 1 4633 E+00 

51 

-.21678E+00 

-•22003E+00 

-•22383E+00 

• 15817E+00 

• 15 113E+00 

•14430E+00 

52 

-• 22255E+00 

— • 22609E +00 

-•23010E+00 

• 15 190E + 00 

• 14696E+00 

• 14 143E+00 

53 

-•22744E+00 

-.23094E+00 

-• 23466E +00 

• 14734 E+00 

• 14320E+00 

• 13850E + 00 


N> 

yo 



UFAC 


I 

JMX 

KMX 

ISUP 

NS 

DPHIMX 

10 

68 

-13 

0 

108 

• 18612E-01 

20 

59 

-10 

0 

109 

•67615E-02 

30 

49 

-9 

0 

123 

•35451E-02 

40 

45 

2 

1 

136 

•22124E-02 

50 

31 

-14 

0 

136 

• 11244E-02 

60 

22 

-15 

0 

136 

•74607E-03 

70 

31 

1 

0 

136 

• 52888 E-03 

80 

31 

1 

0 

136 

• 3922 8E-03 

90 

31 

1 

0 

135 

• 30070E-03 

100 

31 

1 

0 

135 

•22584E-03 

110 

31 

1 

0 

134 

• 16935E-03 

120 

31 

1 

0 

134 

. 12720E-03 

130 

31 

1 

0 

134 

•95640E-04 

140 

31 

1 

0 

134 

•72079E-04 

150 

31 

1 

0 

134 

• 54323E-04 

160 

46 

3 

0 

116 

• 33826E-02 

170 

32 

3 

0 

114 

• 30474 E-03 

180 

22 

3 

0 

112 

•171 89E-03 

190 

28 

3 

0 

112 

• 81272E-04 

200 

37 

8 

1 

112 

•69777E-04 

210 

39 

7 

1 

111 

• 52 5 32E-04 

220 

27 

2 

0 

111 

•30374E-04 

230 

25 

-3 

0 

111 

• 19915E-04 

240 

27 

-6 

0 

111 

• 1 3692 E-04 

250 

28 

-7 

0 

111 

• 99122E-05 

260 

28 

-7 

0 

111 

.731 5 8E-05 

270 

28 

-6 

0 

111 

•54218E-05 

273 

28 

-6 

0 

111 

•49535E-05 


0. 

0. 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 


•98176E+00 

•98176E+00 

•98176E+00 

•98176E+00 

•98136E+00 

•96091E+00 

•98093E+00 

•98103E+00 

• 98 107E + 00 
•98110E+00 

• 98 112E+00 

• 98113E +00 
•98113E+00 


30 



MODEL PERTURBATION (U AND V ARE PERTURBATIONS FROM UINF BUT NORMALIZED BY UREF) 
JP « 44# AMINF - » 76745 > LINF/URfF - . 90113# . CP STAR - -.53345E+00# OFAC • .97372 


X< J-.5) 

YMOD 

C P ( T N F ) 

U 

J 

X(J ) 

V 

DELTA PHI 

DELTA V 

-.00040 

o.oocco 

.80657 

-.42744 

31 

.02400 

.65838 




o.occoo 

1 .03564 

-.67458 



-.60217 

•01206 

1.26054 

.04840 

.03213 

-.68965 

•45648 

32 

.072 80 

.24158 




-.02939 

.11051 

.02079 



-.20711 

.03332 

.44870 

.09720 

.04392 

-1.04515 

.51950 

33 

.12160 

.11919 




-.03949 

-.08506 

.05788 



-.14675 

.05505 

.26594 

.14600 

.04973 

-1.09732 

.52721 

34 

.17040 

.06080 




-.04665 

-.17875 

.09478 



-.10103 

.07695 

.16983 

.19480 

.05309 

-1.09649 

.52212 

• 35 

.21920 

.03744 




-.05158 

-.23178 

.11538 



-.06503 

.09680 

.10247 

.24360 

.05492 

-1.10594 

.52502 

36 

.26800 

.00981 




-.05476 

-.25182 

.12244 



-.03664 

.11645 

.04665 

.29240 

.05540 

-1.11334 

.52624 

37 

.31680 

-.01070 




-.05656 

-.25974 

.12502 



-.01090 

.13612 

.00020 

.34120 

.05468 

-1.11461 

.52914 

38 

.36560 

-.02608 




-.05709 

-.25209 

.12108 



.01471 

.15604 

-.04060 

.39000 

.05360 

-1.12612 

.53553 

39 

.41440 

-.04382 




-.05637 

-.21860 

.10561 



.03706 

.17702 

-.08088 

•43880 

.05146 

-1.15136 

.54950 

40 

.46320 

-.06552 




-.05456 

-.16913 

.08288 



.05473 

.19979 

-.12025 

.48760 

.04827 

-1.17575 

.56392 

41 

.51200 

-.08614 




-.05189 

11920 

.05998 



.07149 

.22438 

-.15763 

.53640 

.04406 

-1.20219 

• .5 802 3 

42 

.560 80 

-.10017 




-.04640 

-.06314 

.03416 



.08936 

.25103 

-.19753 

.58520 

.03879 

-1 .21611 

.59062 

43 

.60960 

-.12456 




-.04404 

.01320 

-.00178 



.10365 

.27994 

-.22821 

.63400 

.03271 

-.76546 

.36253 

44 

.65840 

.03544 




-.03898 

.10576 

-.04708 



.10876 

.29993 

-.07332 

.68280 

.03444 

-.45206 

.21549 

45 

.70720 

-.10174 




-.03368 

.19998 

-.09544 



.10166 

.31510 

-.20340 

.73160 

.02947 

-.34401 

.17282 

46 

..75600 

-.13789 




-.02871 

.28313 

-.14046 



.07990 

.33039 

-.21779 

.78040 

• C22 74 

-.22692 

.12122 

47 

.00480 

-.15458 




-.02482 

.33837 

-.17189 



.05161 

.34469 

-.20619 

*82920 

.01520 

-.12857 

.07604 

48 

.85360 

-.16466 




-.02230 

.37128 

-.19113 



.02267 

.35773 

-.18733 

.87600 

.00716 

-.02881 

.02804 

49 

.90240 

-.16498 




-.02119 

.39884 

-.20697 



-.02059 

.36920 

-.1443.9 

.92680 

-.00089 

.05916 

-.01631 

50 

.95120 

-.15909 




-.02220 

.39004 

-.20113 



-.07600 

.37822 

-.08309 

.97560 

-.00865 

.14676 

-.06261 

51 

1.00000 

-.14183 




-.02590 

.26031 

-.14014 



-.07254 

•36200 

-.06929 

1.02440 

-.01557 









-.02944 



VELOCITIES 

XV 

• 78165E + 01 

• 64734 E + 01 
•54660E+01 

• 4 682 5 E + 0 1 

• 40 5 57 E +01 
•3542FE+C1 
•311 55E+C1 
•2753FE+01 
•24439E+01 

• 21753E +01 
•19402E+01 

• 1 7 32 8 F +01 
•15484E+C1 

• 1 30 3 5 E + 01 
•1225CE+01 

• 11 007 E + 01 

• 97 862 E +00 
•86713E+0C 
•76493F+00 
•67G91E+0C 

• 58 4 1 2 E + 00 
•50376E+CC 

• 4291 A E + 00 

• 35 967 E +00 
•29483E+00 
•23417E+00 
•17730E+00 
•12388E+00 
•73600E-01 
•2490CE-01 
•24000E-01 

•7260QE-01 

•12160E+00 

* 17040E + 00 

•21920E+00 

•26800E+00 

•31680E+00 

•36560E+00 


ON CENTER LINE 
VF 

• 4 83 34 E-02 
•58C19E-02 

• 6 6 1 95 E "02 
•78896E-02 

• 90 129F-C2 
•101Q0F-01 
•11421F-C1 
•12709F-01 

• 1 40 58 E-G 1 

• 15 4 73E-01 
•16960F-01 
•18528 F -01 
•20ie4F-01 
•21941F-C1 
•23811 E-01 

• 25 81 0 F— 0 1 
•2 79 5 5 F-C 1 
•30270F-01 

• 327 84 E-0 1 
•35531 E-01 

• 3 8560E-01 
•41936F-01 

• 4574 9 E-01 
•50134E-C1 

• 5 5298E-0 1 

• 6 1 5 90 E-01 

• 6964 5 E-01 
•P0749E-01 

• 97794 E-01 
•i2eeeE+oo 
•65838E+00 

-•60217E+00 

• 241 5 8F + 00 
-•2071 1E+00 

• 1 1919E + 00 
-•14675E+00 

•68604E-01 
-•10103F +00 
•37442E-01 
- • 6 5 C3CE-01 
•98117E-02 
- • 366 39 E-01 
— • 10704 F-01 
- • 10904 E-01 
-•26003 F-C 1 
• 1 471 5 E-01 


VW 

-•41 861E-02 
— • 465 16 E— 02 
-•48291E-02 

— •5 0390E-02 
-• 53 54 8 E-02 
— • 5 76 0 4E-02 
-•62175 E-02 
-• 66 894 E-02 
-•7149 3 E— 0 2 
-•7 5 8 10 E-02 
-•79768E-02 
-• 83342E-02 
-• 8 6544 E-02 
-• 8 93 99E-02 
-•91941E-02 
— * 94 203 E-02 
-• 962 17E-02 
-• 9 80 13E — 02 
— • 9961 5 E— 02 
—•10 1 04 E-01 
-•10232 E-01 
— • 10346E-01 
-• 1 04 47 E— 01 
-•10537E-01 
-• 1061 8 E-01 
-•10 6 89 E-01 
-• 10751 E-01 
-•10 805 E-01 
— • 108 52 E— 01 

— • 10902 F — 01 
-•113 94 E— 01 

— •11245E— 01 
-•11 1 05F-01 

— • 11043 E-01 
-•11 303E-01 

— • 11245E-01 
-•113 50E-01 

— •112 94 E-01 
-•11 3 71 E-01 

—•113 17E-01 
-• 11402E-01 

-.11349E-01 
— • 11 4 29 E-01 

-.11376E-01 
-• 11448 E— 01 

— • 11395E-01 


XU 


-•71489E+01 

— • 59 719 E + 01 
-• 50756E+01 
-. 43700E+01 

— *379988 + 01 
-• 33296E+01 
-.29349E+01 

— • 25 991E + 01 
-• 23097E+01 

— •20578E + 01 
- • 10 366 E + 01 
-• 16407E + 01 
-• 14660E+01 
-• 13093 E + 01 
— • 11679E+01 
-• 10397E+01 

— • 922 90 E + 00 
-• 81 606 E +00 

— • 71794E + 00 
-.62753E+00 
-. 5439 6 E +00 

— •46647 E + 00 
-.39442E+00 
-* 32726E + 00 
-• 26451 E+00 

— • 20 5 7 4 E + 00 
- • 1506 OE + 00 
- • 98 745 E-01 

— • 49200 £-0 1 
- • 40000 E-03 

• 48 400E-01 
• 97200E-01 

• 14600E + 00 

• 194 80E + 00 

• 24 3 60 E + 00 
• 29240E+00 

• 341 20E + 00 


UF 


-.10242E-02 

— •13779E-02 

— • 18719 E-02 
-.24890E-02 

— •32202 E-02 

— • 40606 E-02 
-.50099E-02 
-.60725E-02 

— • 725 80E-02 

— • 85807E-02 

— *10060 E-01 
-•1 1722E-01 
-.13596E-01 

— • 15723 E-01 
—•181 51 E-01 

— • 20939E-01 
-.24164E-01 

— • 279 23F-01 

— •32343 E-01 

— • 37594E-01 

— • 43908 E-01 
— • 51 611 E-01 
-•61178E-01 
-.73342E-01 

— • 89294E-01 

— • 11117E + 00 

— •14332 E + 00 
—•1963 1E+00 
-•29954E+00 
-.42744E+00 

— *674588+00 

• 45648E + 00 

• 20793 E-01 
•51950E+00 

•5788 3E -01 
•52721 E+00 

• 94778 E-01 

• 5221 2 E + 00 

• 1 1 5 38 E + 00 
•52502E+00 

• 12244 E +00 
• 52824E+00 

• 1 2 502 E+00 
•52914E+00 

.12108 E+00 


UW 


— *392908-03 
—•18301 E-02 
-•38246E-02 

— • 568 4 9E-02 

— •73341E-02 

— • 8 7452 E-0 2 
-•99330E-02 
-• 1 0932 E-01 
-.11781E-01 

— • 1251 3E-01 
— • 1 3152E-01 
-• 13718E-01 
-•14221E-01 
-• 1 4671 E-01 

15074E-01 
-•1 5435E-01 

— • 15760E-0 1 
-• 16052E-01 
-• 16315E-01 
-• 165 52E-01 

— •16766 E-01 
-• 16959E-01 
-•17130E-01 

— • 17280E-01 
-•1741 1 E-01 
— • 17523E-01 
-•1761 3E-01 
— • 17672E-01 
— • 17637E-01 
-•17131 E-01 

— * 171 31E— 01 
-•1 8290E-01 
— • 18290E— 01 
— • 18435E-01 
— • 1 8435E— 01 
-• 18358E-01 
-•18 3 58 E— 01 
-• 18334E-01 
-•1 8334E— 01 
-•183 72E-01 
-.18372E-01 

— •1842 IE-01 

-. 10421E-O1 
-*18481 E-0 1 
-. 18481E-01 



** * 


.4 1 4 4GE +00 
• 46320E+00 

• 51 200E + 00 

• 5 600 OE + OO 
•60960E+00 
.65840E+00 
• 70720E+00 

• 75600E + 00 
•80480E+00 
• 85360E+00 

• 90240 E +00 

.95120E+00 

.1GC00E+01 

•1048FE+C1 
•10976E+01 
•11479E+C1 
. 12013E+C1 
•12582E+01 
•13196E+01 
•13837E+G1 

• 1453 1 E + C1 

• 1 5 278 E +01 
.160818+01 

• 1 6949 E + 01 
•178P9E+01 

• 1891 1 E + 01 

• 2002 6E +01 
•21247E+01 

• 22 590 E +01 

• 24075 E + 01 

• 25724 E + 01 
•27566E+01 

• 296^ 2E +01 

• 31 993E+01 
•34679E+01 

• 37778 E + 01 
•41395E+01 


-.43P24E-C1 
•37057E-01 
-.65517F-01 
.547358-01 
-.86136F-01 
•714 9 5F -01 
-.106178+00 
.893568-01 
*•124 5 6F +00 

• 1036 5 E +00 
•35439 E -01 

• 1 0676 E + 00 
-•10174E+00 

• 1 0166 r + 00 

— • 1 37 8 9E +00 

• 79 89 9 F-01 

— • 1 54 5 6 F +00 

• 5 1 6 1 2 F-01 
-.16466F+00 

• 22 6 66 F-01 
-.16498E+00 

- • 2C 5 92 F-01 
-.15909r+00 

- .76002 F-01 
-•14163E+00 

- « 72 540 E-01 

— *776 20 8-01 

— • 6 6799F -0 1 
-• 59 28 4 F-01 
-.53443F-01 
-.46652E-01 
-.44579E-C1 
-.41029r-01 
-•37877C-01 

— • 3 5 04CF-01 
-•3 24 59F-01 

— • 30092 E -0 1 

— • 27 90 5 F-01 

— .25871 E--01 
- • 2 3982 F-01 
-.2 220FF-C1 
- • 2 05 43E-01 

— • 1 8974 F-01 

— • 1 7492E-C1 
—•16091 F-01 
-.14761E-01 
— • 1 3498F-01 
—•122946-01 

— • 1114CE-01 
-•10024E-01 


LO 

U> 


39000 E+00 


•1147 6 E-01 

. 39000 E+00 

-. 11 4 22 E-01 

•11 522E-01 

.43880E+00 

- • 1 1 467E-01 


.115798-01 

.48760E+00 

115 22E-01 

• 11 653 F-01 

. 53640E+00 

-.11595 E-01 

• 1 1 745 E-01 

.585206+00 

-.116 85 E-01 

• 118 38 E-01 

• 63400E+00 

-.117768-01 

• 12 660 E-01 

.6 8280E + 00 

-.12553 E-01 

• 1 2 706 C-01 

•73160E+00 

-.12 5 96 E-0 1 


• 12 563 E-01 

. 78040 E+00 

-. 1 2450 E-01 

-.124058-01 

. 82920E+00 

-.122918-01 

■•12281 E-01 

. 878 OOE+OO 

-• 12166E-01 

• 1 21 82 E-01 

• 92 6 8 OE + 00 

-• 12070E-01 

. 12 076 E— 01 

• 97560E + 00 

-.11 968 E-01 

. 11 929 E-01 

•10244E+01 

.118478-01 

• 107 3 2 E + 01 

.117678-01 

• 11227E+01 

•11 685E-01 ■ 

. 11746E+01 

• 1 1 5 99 E-0 1 

• 1 2 297 E +01 

♦ 1 1 5 06 E-01 

•12885E+01 

• 11 404 E-01 

• 1 351 2 E + 01 

. 11 290E-01 

•14184E+01 

.111608-01 

. 14904 E + 01 

. 1 10 12E-01 

• 1567 9E+01 

• 1 0840 E-01 

• 165 15E + 01 

• 10641 E-01 

. 174 1 9 E +01 

• 1 0409 E-01 

.1 8400E+01 

• 101 3 8E-01 

. 19 468 E + 01 

• 96 209 E-02 

• 20636E+01 

• 94 506E-02 

. 21918 E + 01 

• 90 180 E-02 

. 23332E+01 

• 85 169E-02 

• 24 899E + 01 

• 793 71E-02 

• 2664 5E +01 

. 72735E-02 

. 2 8604E + 01 

• 65246E-02 

. 3081 6E+01 

•56978E-02 

.333346+01 

•^81 63E-02 

. 36227E+01 

• 39319E— 02 

.395846+01 


•53553E+00 -.18562E-01 

• 10561 E +00 — • 185 62E-01 

• 5495 0E + 00 — *186458-01 

.82 880 E-01 - • 18645E-01 

• 56392E + 00. -•18715E-01 

• 5 99 76 E-01 — .1 87 15 E-0 1 

•580 23 E+00 18778E-01 

• 34176E-01 - • 1 8778E-01 

• 59062 E + 00 — • 1 8817E-01 

- • 1 7 804 E-02 - • 1 8 8 1 7 E-01 

• 362 53 E +00 -.18861E-01 

— • 4 7076 E-01 -•10 061 E-01 

•21549E+00 -• 19122E-01 

-.95443 E-01 - • 19122E-01 

• 17282F + 00 - • 1 82 94 E-0 1 

-.14046 E+00 -•18294 E-01 

• 12 122E + 00 -• 1 773 9E-01 

-.17189E+00 -• 17739E-01 

• 76041E-01 -• 17442E-01 

-.191 13 E+00 -.1744 2 E-01 

• 2 0O39E-O1 -• 17278E-01 

-•20697E+00 -.17278E-01 

-•16313 E-01 -.17186 E-01 

— • 20113E + 00 -.171 66 E-01 

- ,626098-01 — • 1 7144E-0 1 

-.1401 4 E+00 — • 17144E-01 

-•111 35 E+00 -• 1 7100E-01 

-.765548-01 -. 17002E-01 

-• 59943 E-01 -.16911 E-01 

-.49094 E-01 -.16 812 E-01 

-.41081 E-01 -.16701 E-01 

-.34798 E-01 — * 16574E-01 

— • 29707E-01 - . 1 6427E-01 

-• 2 5494E-01 -.1625 8 E-01 

-.219 5 7E-01 -.16062 E-01 

-.189 56 E-01 -• 1 5836E-01 

-.163 91 E-01 -.15 574 E-01 

- . 14 18 5 E-01 -. 1 5273 E-01 

-• 1 22Q0E-01 -• 14925 E-01 

— . 1 0630E-01 -.145268-01 

— . 91 96 9E-02 -.14066 8-01 

- .795 13E-02 -• 1 3536E-01 

— .6 8686E-02 -.12 924 E-01 

-.59292 E-02 -.12218E-01 

- • 5 1 172E-02 —.11 3 99E— 0 1 

-•44211E-02 -• 104 44E-01 

— • 38337E— 02 -.9323 5 E-02 

-.33531 E-02 -. 79950E-02 

— .2 9054E-O2 -.6398 7E-02 

- . 27497E-02 -.44478E-02 



I 

J MX 

KMX 

IS UP 

NS 

DPHIMX 

DA 

X 

52 

'1 

0 

111 

-.387705-02 

0. 

2 

52 

-1 

0 

111 

* 37 16 38-02 

0. 

3 

52 

-1 

0 

111 

-•30487E-O2 

0. 

4 

52 

-1 

0 

111 

•266405-02 

0. 

5 

52 

-1 

0 

111 

-•211 905-02 

0. 

6 

5? 

*1 

0 

111 

• 178055-02 

0. 

7 

52 

-1 

0 

111 

-.13820E-02 

0. 

6 

53 

-1 

0 

111 
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